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We investigate the large scale anomalies in the angular distribution of the cosmic microwave 
background radiation as measured by WMAP using several tests. These tests, based on the multipole 
vector expansion, measure correlations between the phases of the multipoles as expressed by the 
directions of the multipole vectors and their associated normal planes. We have computed the 
probability distribution functions for 46 such tests, for the multipoles £ = 2 — 5. We confirm earlier 
findings that point to a high level of alignment between 1 = 2 (quadrupole) and £ — 3 (octopole), but 
with our tests we do not find significant planarity in the octopole. In addition, we have found other 
^~>' possible anomalies in the alignment between the octopole and the £ — 4 (hexadecupole) components, 

2 ', a-s well as in the planarity oi £ — 4 and £ — 5. We introduce the notion of a total likelihood to 

estimate the relevance of the low- multipoles tests of non-gaussianity. We show that, as a result of 
these tests, the CMB maps which are most widely used for cosmological analysis lie within the '^ 
10% of randomly generated maps with lowest likelihoods. 

PACS numbers: 98.80.-k, 98.65.Dx, 98.70.Vc, 98.80.Es 
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rn ■ I. INTRODUCTION 

o , 

\^ , The cosmic microwave background (CMB) anisotropies have been measured with exquisite accuracy by WMAP, 

f^ . and the impact on Cosmology has been profound [IQjyiH- However, as the ACDM cosmological model becomes 

standard lore and the parameter space becomes narrower, the focus naturally drifts to the apparent anomalies. Among 

sources of concern that have survived the WMAP 3-year data are the lack of large-angle correlations [H Q, |j] , which 

is mainly due to the low value of the cosmic quadrupole |j, la, S S S] j and the alignment between the quadrupole 

O ; {i = 2) and octopole {£ = 3) He, 9, 10, 11, 12, 13, 14, ISJ. 

-*— > The combined statistics of these effects, which were already in the COBE data [iajll3, implies that the probability 

jrt , that our CMB sky was generated in a random process is only 0.005 - 0.02%, depending on the map that is being 

!• ■ tested. 

. ^ I Several recent works have reported some anomalies in the data: the low-order multipole values [E Q, II, H, El, El, 

b>^ . El, llll ; the alignment of some low-order multipol es la. llSl l2ll 122 . l23 | ; an unexpected asymmetric distribution on 

Vh ' the sky of the lar ge-s cale power of CMB data [Ij, lU |2l |2l |23| ; indications for a preferred direction of maximum 

5t , asvmmetrv |^ . |2 5l. I27l l2a . 29, ^(J|; as well as apparent non-gaussian features detected via the wavelet method or other 

analyses p3ll 133 . 133 . l34. 35 . 36] . These anomalies have motivated many explanations, such as compact topologies 

|37l l38| , a broken or suppressed spectrum at large scales |3l HI llj, IH |4l |4J| , oscillations superimposed on the 

primordial spectrum of density fluctuations J45l l4a . |4 ^. anisotropic cosmological models |4l, |4l, |5l| and possible 

extended foregrounds that could be affecting the CMB [10, 51, 52, 53, 54] . 

In this article we examine the multipoles £ — 2 — 5 and search for anomalies in their phase correlations. We measure 
these correlations through alignments between the multipole vectors or through their associated normal planes |9(. 
We also look for evidence of planarity (or self- alignments) in each individual multipole, and for evidence of alignments 
between the multipole and normal vectors with some specific directions in the sky, such as the dipole axis, the ecliptic 
axis and the Galactic poles axis. In total we have considered 38 tests of alignments between multipoles, plus 8 tests 



Oh' 



•Electronic address: abranioQfma. if. usp. br| 
^Electronic address: bernui@das.inpe .br| 
-f Electronic address: ivanOdas.inpe.bn 
§ Electronic address: thyrso@das.inpe br| 
^Electronic address: alex@das.inpe.br| 



of alignments of the multipoles with a priori directions. We have computed the probabihty distribution functions 
(PDFs) for those tests using 300,000 mock maps. 

For the statistical analyses peformed here, we need the aim's, (£ — 2 — 5) of each CMB map under investigation. The 
corresponding a^m's were extracted, after applying the Kp2 WMAP mask (to minimize foreground contaminations 
from the beginning) , through the HEALPix routines ,5^ . Our statistical tools were then used to analyze the WMAP 
1-year and 3-year data Internal Linear Combination maps 0, 3 (henceforth ILC), the co-added 1-year and 3-year 
WMAP data, as well as the cleaned CMB maps of Tegmark et al. jlSl Il5l | based on 1-year WMAP data (henceforth 
TOH). 

This paper is organized as follows. In Section II we summarize the multipole vector formalism and the several 
different statistics that can be used to test for alignments and phase correlations within a given multipole. In Section 
III we briefly describe the CMB maps used. Section IV presents our statistical tests and the results of the PDF 
computations for those tests. We also analyze the salient features of the CMB maps and discuss which tests can be 
considered anomalous, and of those, which are robust and which are most sensitive to noise. The conclusions are 
presented in Section V. 

II. MULTIPOLE VECTORS, NORMAL VECTORS AND STATISTICS OF PHASE CORRELATIONS 

Multipole vectors were introduced in CMB data analysis by Copi et al. 9] , and Katz and Weeks [l^, ll^ found 
an elegant algebraic method to compute the multipole vectors given the spherical harmonic components aim ~ see 
also iSw] for an alternative algebraic method and |6l ll0lll4| for purely numerical methods. The multipole vectors are 
essentially eigenvectors - i.e., they are solutions of a set of polinomial equations whose parameters are the multipole 
components agm- 

The idea of the multipole vector representation goes back to J. C. Maxwell in the XIX**^ century: the multipole 
decomposition of a field f{9, 4>) on S'^ implies that for each multipole (. there are (. eigenvectors of norm unity, n^^'P' . 
Since there are only 2£ phases for each multipole, the spherical harmonic representation and the multipole vector 
representation have the same number of degrees of freedom in each individual multipole: 

^lip^^ ^ aimYemi9,v) ^ Dil[n('^P^ .h{0,^) - Zi_,{e,ip) , (1) 

m——£ p—1 

where Zi-i just subtracts the residual £' < £ total angular momentum parts of the product expansion, and is irrelevant 
to our analysis - see |ll| for an enhanced discussion of the multipole vector expansion. 

It can be seen from the product expansion above that, whenever using the multipole vectors to test for alignments, 
it is irrelevant what the amplitudes of the multipoles are - just their phases matter. This is the main feature of the 
tests based on the multipole vectors that sets them apart from other tests of non-gaussianity. 

Notice that, contrary to the C^'s, which are always positive-definite, the D^'s of Eq. (Q can be either negative or 
positive. Because of the product expansion in the right-hand-side of Eq. |^, switching the sign of De is equivalent to 
switching the signs of an odd number of multipole vectors, and switching the signs of an even number of multipole 
vectors leaves the sign of Di invariant. Therefore, the product expansion in Eq. 1^ has a sign degeneracy in the 
amplitudes Dg as well as in the the multipole vectors n^^'P\ 

This "reflection symmetry" n^^'^-' <-> — n^^'^' implies that the multipole vectors define only directions [Ilj, hence 
they are "vectors without arrowheads" living on the half-sphere with antipodal points identified, or S^/Z2. This space 
is also known in the literature as the real projective space MP^, and is useful in the characterization of nematic liquid 
crystals, where the orientation of the molecules is an order parameter - but it makes no difference where heads and 
tails are jsj- 

It is extremely useful to represent these directions as vectors, but for that we will need to cope with the degeneracy 
in representing these directions. For each multipole order £ there is a 2^~^-fold degeneracy in the signs (or orientations) 
of the multipole vectors - corresponding to the £ signs of the multipole vectors that can be switched arbitrarily, divided 
by two to account for an irrelevant overall sign which is determined by Di. We can break this degeneracy by always 
working in one particular hemisphere, and any such choice will automatically determine the signs (orientations) of all 
multipole vectors - as well as the sign of the D^'s. However, we should be aware that in doing so we are necessarily 
picking one of the 2^~^ possible sign conventions in the product expansion of Eq. Q. As we will discuss below, this 
is not a problem as long as we use invariant tools which are not sensitive to the signs of each multipole vector. 

Starting with the £ multipole vectors one can also construct £{£ — l)/2 = A normal vectors - or normal planes - 
defined as: 
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•^■« = h'^'P A h'^'P' , (p / p' , g = 1 . . . A) . (2) 



By construction, because the multipole vectors define only directions, the normal vectors also possess reflection 
symmetry, w^^^i' ^^ —w^^'ii . But the normal vectors need not be (and generally are not) of norm unity, so instead 
of living in S"^ jTL^ the normal vectors belong to the space ]R'^/Z2 - which is isomorphic to S0(3), see [S3- We can 
still break the degeneracies imposed by reflection symmetry by defining all normal vectors so they lie in the same 
hemisphere as the normal vectors. However, just as before, we should be aware that in doing so we are choosing one 
of many possible representations for the normal vectors. 

To summarize, for £ = 2 there are 2 multipole vectors, n^^'^' and n*-^'^), and only one normal vector, w^'^'^' = 
fj(2,i) X fi(2,2). £qj. £ _ 3 there are 3 multipole vectors and 3 normal vectors; and so forth. These constitute the basis 
for the statistical tests defined below. 

A. Properties of the tests under reflection symmetry 

We will define below, in Sec. IIC, a series of tests which are manifestly invariant under the refiection symmetry 
n <-^ —h which characterizes RP^ (where the multipole vectors live) as well as R^/Z2 (where the associated normal 
vectors live.) 

Our motivation for this remark is that if a given test is not invariant then its validity and usefulness is questionable. 
In particular, one should be careful not to employ tests which depend on the choice of hemisphere to represent the 
vectors. This seems to be the case of some of the tests that have been used to estimate the "planarity" of the CMB 
maps, if these tests make use, in one way or another, of the notion of "average vectors" . The reason there is no such 
thing as an "average multipole vector" or an "average normal vector" is simple: the "vector sum" operation does not 
yield a singly valued result due to the refiection symmetry. In fact, the result of "summing" two vectors of ]R'^/Z2 
would be a degenerated pair of directions: 

(±.i,»(±« = {±:;;±* . (3) 

In general, by "summing" fc vectors of M'^/Z2 one obtains 2*^"^ vectors, corresponding to the 2^ possible permutations 
of the ± signs of each vector, divided by two to account for the symmetry v •;-> —v of the resulting vector. This means, 
in particular, that the "average multipole vector" is in fact an object 2^~^-times degenerated, and that the "average 
normal vector" is an object 2^(^~^'/^^^-times degenerated. 

Obviously, by fixing a hemisphere to represent all vectors one breaks this degeneracy, but this just hides the plain 
fact that by doing so one is simply choosing (rather arbitrarily) one of many possible representations, and one of 
many possible answers for the sums of those vectors. Consequently, unless these degeneracies are properly taken into 
account (by, e.g., symmetrizing over all possibilities or ordering the results by norm), any test which is derived from 
the notion of summing multipole or normal vectors is fiawed. 

B. Global estimates of non-gaussianity 

We will investigate large-scale correlations in the CMB maps within single multipoles and between different multi- 
poles by measuring the alignments between either the multipole vector themselves, or between their associated normal 
planes. To be sure, there is no upper limit to the number of tests we can devise to search for non-gaussianities, and in 
testing any fixed sample such as the CMB one should bear in mind that there are always some statistical tools which 
will yield a positive detection given some arbitrary criteria. 

In practical terms this means that if we perform a large number TV of independent tests on a fixed sample, then 
we should treat the results of these tests themselves as random numbers. Therefore, when performing many tests on 
a map and searching for clues of non-gaussianity one must always look at the complete set of results for those tests 
and at the total probability (or some estimate of the likelihood) that map is a realization of a random process. If 
that likelihood turns out to be very small compared to the typical likelihoods of gaussian maps, then one can look 
for the particular test (or tests) that is likely responsible for that anomaly. If, however, a certain test turns out to be 
"suspicious" but the likelihood is not anomalously small, then we cannot rule out the possibility that the result for 
that particular test was just a fluke. 

Assuming that random processes are indeed behind the mechanism that generated the sample, we can define a total 
likelihood in the following sense. Suppose we have TV statistical tests such that each test Ti (i = \ . . .TV) is a random 
number in the interval < T^ < 1, with normalized probability distribution functions Pi(Ti). Given a sample (a map 



M) with Ti — T^ , the probabihty that a random sample has a value of Ti higher than T^ is 



P,+ {Tl') = f dTP,{T) , (4) 



and the probability that a random sample has a value of Ti lower than T^^ is: 

n-iTt") = / dTP,{T) = 1 - P+iTt') . (5) 

Jo 

Evidently, for the median value Ti we have Pi+{Ti) = Pi-{Ti) = 1/2. We will define the likelihood of the map M, 
given the TV tests, to be: 

L^M) = Y[2P+{T,^) X 2P._(r^) = A^l[P.+ {Tt') [l - P.+ m^)] , (6) 

1=1 i=l 

where the factors of 2 have been inserted for normalization purposes, in order to make L^{M) = 1 for a map M 
whose tests are all exactly equal to their median values. This likelihood measures the total probability that the map 
M does not have too low and too high values of the tests Ti. Evidently, any deviation of the tests from the medians 
will decrease Lj^/. 

Of course, Eq. © is itself an arbitrary definition, and we might as well have used the expectation values instead 
of the medians to define the likelihood. Indeed, one could switch the roles of the medians with the expectation values 
in the procedure above and still our results would be very similar. 

It should be stressed that the total likelihood defined in the sense above should not be interpreted as the probability 
that a particular map was generated by a gaussian mechanism - it is merely an estimator of how much that particular 
map deviates from a typical one, given the J\f statistical tests of non-gaussianity. In Sec. IID we construct such 
likelihoods, and compute their distributions assuming random phases. 

C. Statistical tests of isotropy 

We now define the statistical tests which will be employed in our analysis of CMB data. We have ensured that all 
tests are invariant under reflection symmetry, so it makes no difference which hemisphere one chooses to represent 
the multipole and normal vectors. 

The tests have been normalized so that they always fall in the interval < Ti < 1. We have generated 3 x 10^ 
simulated (mock) maps using gaussian random phases, and the resulting PDF's for the tests are shown in Figs. 1-5. 

1. S statistic 

The S statistic is a widely used tool [l^ [lj| , and it measures the alignment between normal planes of different 
multipoles. It is defined as: 

q=lq' = l 

where, as defined above, A — £{i — l)/2. We can also use S in just one multipole, in which case the normalization is 
a bit different: 

5„.-^x: K-^-^'^-'t (8) 

q,q'>q 

The statistic Su measures the "self-alignment" of the normal vectors, and is related to the "planarity" tests p.[l^. 
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FIG. 1: Normalized PDF's for the S statistic found by simulating 3 x 10^ mock maps, binned in intervals of 0.01. From left 
to right, top to bottom: S22, S23, S24,, S25, S33, S34,, S35, S44,, S45 and Sss- In all panels, the horizontal axis correspond to the 
value of each individual test, and the vertical axis to its normalized PDF. 



2. D statistic 
This is analogous to the S statistic, but the D test disregards the norm of the normal vectors: 
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q=\q' = \ 



^i' 



(9) 



We can also use the D statistic within a single multipole, as was done for S. However, this test only gives nontrivial 
information for ^ > 3. With the proper normalization we have: 
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(10) 
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FIG. 2: Normalized PDF's for the D statistic. From left to right, top to bottom: D23, D24, D25, D33, D34, D35, D44, D45 and 
D55. 

3. R statistic 



A similar tool is the R statistic, which measures alignments in essentially the same way as the S statistic, but it 
uses the multipole vectors instead of the normal vectors: 



nii,P) . f,(i\p') 



p—i p'—i 
Within a single multipole, the R statistic is suitably defined as 
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(11) 



(12) 



4- B statistic 



We can also test if the multipole vectors align with the normal vectors. Hence we define the B statistic: 



^"'-iEEl-*'^'^-^"'^'^ 



, e^i' 



q=lp' = l 

Within a single multipole, the B statistic only gives nontrivial information for ^ > 3, and we have: 
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(13) 



(14) 



Notice that, as opposed to the 5*, D and R statistics, the B statistic for I ^ i' \s, not symmetric, Bu' 7^ 5^'^- For 
simplicity, in the present approach we have only considered the cases B^i where I < t' . 
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FIG. 3: Normalized PDF's for the R statistic. From left to right, top to bottom: R22, ^23, R24, R25, R33, R34, R35, -R44, Rib 
and R55. 

5. N statistic 
We can test if the multipole vectors align in a particular direction Z by using the N statistic: 



N,^'-±\ni^^^^.Z 



0=1 



(15) 



6. W statistic 
We can also test if the normal vectors align in a particular direction Z by using the W statistic: 



vt/,^i^L(^^^).z 



q=l 



(16) 



Note that the test S'-^^'^^ of Ref. Hi is a combination of W2 and W3, namely: S'-^-'^^ = {W2 + 3W^3)/4. 
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FIG. 4: Normalized PDF's for the B statistic. From left to right, top to bottom; B23, B24, B25, -B33, B34, B35, B44, B45 and 

B55- 
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FIG. 5: Normalized PDF's for the A'' and W statistics - alignments of the multipole (A'') and normal (W) vectors with a 
particular direction in the sky. Top line, from left to right: A^2, N3, N4, N5. Bottom line, from left to right: W2, W3, W4, W5. 



D. Likelihoods 



With the tests defined above we can now compute likelihoods as in Eq. ©. However, we have decided not to include 
the alignment tests N and W in the analysis, as they test correlations with an a priori direction, which we find rather 
arbitrary compared to the other tests. Nevertheless, it should be noted that there are significant correlations between 
the quadrupole and the octopole with both the ecliptic plane and the direction of the cosmic dipole. Although these 
correlations appear to be too strong or too weak only at > 95% C.L. as measured by our tests Ni and W\, when 
combined in the statistic S^^'^^ of Copi et al. [23 ■ we obtain a result which is > 99.5% C.L. These correlations have 
been treated in much greater detail in Refs. |il Il4l l2lll27l | . 

Using the 38 tests S, D, R and B above, we can define the total likelihood of a map, L^siS, D, R, B), as in Eq. 10. 
We have computed L^g for 25,000 random maps, and we show a normalized histogram for Logio.^38 in Fig. 6 (left 
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FIG. 6: Normalized histograms of Logio 1/38(5', D, R,B) (left panel), LogioL29(5, D, R) (center) and LogioI/29(5', i?, B) (right), 
obtained by simulating 25,000 mock maps. 

panel.) 

Since we have 38 tests but only 28 independent random degrees of freedom in the i — 2 — 5 muhipoles, we have 
found useful to define hkelihoods using a subset of the complete set of tests. We have thus defined the likelihoods 
L29{S, D, R) using the tests S, D and R, and ^29(6', R, B), using the S, R and B tests. Their normalized histograms 
are also shown in Fig. 6 (center and right panels.) 

III. CMB MAPS 

In this work, we use five WMAP CMB maps (three derived from the 1-year data [l| and two derived from the 
3- year data |53|): the 1-year and 3-year co- added maps [53|; the 1-year and 3-year ILC maps 0,|5a|; and the TOH 
map _l5j. 

The co-added WMAP map results from the combination of the eight differential assemblies (DA) [53 in the Q-, V-, 
and W-bands, listed above, using the following inverse-variance noise weights method. Thus, for any co-added map, 
the temperature of pixel n is given by: 



Tin) 



Ell T.in)/<Jnn) 



where Ti{n) is the sky map for the DA i with the foreground galactic signal subtracted, and where: 



(17) 



(18) 



is the pixel- noise per observation for DA i. The eight values of ctq ^ are the pixel-noise for DA i, and can be found for 
each set of released maps in ^(l| . 

The WMAP 1-year Internal Linear Combination map (ILC-lyr) Q has been built mostly to convey a visual 
information of CMB anisotropics. It is composed by the five WMAP temperature intensity maps, through a weighted 
linear combination, to minimize Galactic foreground contamination in twelve regions of the sky, eleven of which lie 
within the Galactic Plane and one which lies outside it. The ILC-lyr map is not reliable for quantitative CMB 
analysis, but in this work we use it only for the completeness of our analysis. The 3-year WMAP ILC map (ILC-3yr) 
158], however, brings some improvements over the ILC-lyr, mainly by dealing with the regions selected for Galactic 
foreground estimates, which make it a reliable estimator of the CMB signal for large angular scales (> 10°) and, 
therefore, suitable for our analysis since we are interested in ^ = 2 — 5. 

The TOH map il^] was constructed with no assumption about the CMB power spectrum, foregrounds or noise 
properties. The only consideration was the Planckian nature of the CMB temperature spectrum. This map is formed 
by the combination of the five WMAP bands considering weights that depend both on angular scale and distance to 
the Galactic plane. The cleaning process is done in multipolc space a;™, weighted also by the beam function of each 
channel (W-band has four channels, Q- and V-band have two channels each, and K- and Ka-band have one channel 
each). They obtain weighting coefficients similar to those used by the WMAP team on large scales. 



IV. STATISTICS OF LARGE-ANGLE ANISOTROPIES 



Given a CMB map, the harmonic components aim can be extracted using HEALPix |5q |. and the multipole vectors 
and their statistics can be easily computed. In this analysis we will use the five maps described above, plus the 
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TOH map cleaned with the mask "M6" given in ^^ [13 ■ For all maps we have used the Kp2 mask based on 3-year 
WMAP data l^i^], then we remove their residual monopole and dipole components. It should also be noted that the 
relativistic Doppler correction to the quadrupole is an important factor that must be subtracted from the maps, since 
it corresponds to a non-primary source of the quadrupole |l4l l27l | . 

In Table I we present an abridged version of the tests that were defined in Sec. IIC. For an easier interpretation of 
the results, we show the probabilities P+{Ti) that a random map would have higher values for the test Ti than the 
value for that test in the actual CMB maps. So, for instance, for the test S22 and for the Co-added 1-yr map we 
quote the value P+(S'22) = 0.294. This means that the probability that a random map has a value of S22 which is 
higher than the value obtained for the map Co-added 1-yr is 29.4%. Therefore, in Table I, any values which are too 
close to either zero or one should be viewed with suspicion. 

It can be immediately seen from Table I that the quadru pole -octopole alignment (revealed by iS'23, I?23 and ^23) 
is robust in all maps, as has been noted by p. l9l lla ll4L l20l l2ll l27l| . Our results confirm that the probability that a 
random map has higher quadrupole-octopole alignment, as measured by £'23, is approximately 1-2%. 

Table I also reveals (through the self-alignments £'33, I?33, i?33 and S33) that, at least for our set of invariant tests, 
there is no evidence that the octopole is significantly "planar" , in the sense defined by these tests. 

There are other anomalies for £ > 3 as well: we find that the octopole and hexadecupole (£ = 4) are mis-aligned 
to a very significant degree, with a probability in the range 1-4% for the test 6*34 and 5-20% for the tests I?34 and 
i?34. The triandabipole {(. = 5) also seems substantially mis-aligned with the octopole i. ~ i, with a probability in 
the range 4-8% for the test S^^. These results, plus the fact that the o ctop ole is significantly aligned with the dipole 
axis, give further support to the conjecture known as the 'axis of evil' |21| . 

Another apparent anomaly that is revealed by Table I is the self-alignment of the multipoles £ = 4 and £ = 5, 
indicated by the values of P+ for 5*44, D55, i?44 and R55. 

Not shown in Table I are the results for the tests N and W, using as a priori directions to be tested against the 
dipole axis, the ecliptic plane and the galactic plane. The results for these 24 tests indicate that there are si gnif icant 
alignments of the quadrupole and octopole with the direction defined by the cosmic dipole 0, I3, [13, 111, 120, 12^ U^ ■ 
W2 ^^° '^ ~ 0.10 — 0.01 and W^ ^^° ° ~ 0.03 — 0.05 depending on the choice of map. We also confirm the results of Copi 
et al. H, llJ, llfl I23 , who found strong correlations of the quadrupole and octopole with the ecliptic plane: we get 
^Ecliptic ^ Q g3 _ g gg ^^^ ^Dipole _ q gg _ g gg jepeuding ou thc cholce of map. If combined into a single test, the 
tests W2 and W3 give rise to very significant correlations, both with the cosmic dipole and with the ecliptic plane. 
For all other multipoles and directions we have found no significant alignments using our tests. 

We can estimate the total likelihoods for the CMB maps of Table I. This is shown in Table II for the likelihoods 
L3s{S, D, R, B), L29{S, D, R) and L2g{S, R, B). Comparing the values of Table II with the normalized histograms of 
Fig. 6, we can see that, apart from the TOH map with the Kp2 mask applied, which lies within the ~ 30% of random 
maps with lowest likelihoods, all remaining maps fall within the ~ 10% of random maps with lowest likelihoods. The 
likelihoods are particularly small for the 1-year ILC map, which, as discussed, is probably contaminated by residual 
foregrounds. 

V. CONCLUSIONS 

We have applied 38 tests of non-gaussianity, as well as 8 tests of alignments with 3 distinct preferred directions, on 
the most widely used CMB maps based on 1-year and 3-year WMAP data. In order to properly analyze the set of 
tests performed in the WMAP maps, we have introduced the notion of a total likelihood to estimate the relevance of 
the low-multipoles tests of non-gaussianity for each map. A cautionary note is sounded by the total likelihoods: this 
criterium shows that the CMB maps we studied have rather low likelihoods, but they still lie within the '^ 10 % of 
random maps with lowest likelihoods, meaning that 90% of the gaussian mock maps have higher likelihoods than the 
CMB maps. 

We confirm the well-known alignment between the quadrupole and the octopole, and we found other significant 
levels of alignment between the octopole and the hexadecupole. Moreover, we detected that the hexadecupole and 
the (. — b multipole have significantly low levels of self-alignments (see Table I) . Taken individually, these tests show 
non-gaussian features at more than 95% C.L. On the other hand, no significant evidence of self-alignment (which is 
related to the planarity) was found for the octopole. 

Regarding correlations with the directions of the cosmic dipole, the ecliptic plane or the galactic plane, we have 
confirmed the correlations of the quadrupole and octopole with the cosmic dipole, as well as the correlation of the 
octopole with the ecliptic plane. For all other multipoles and directions we have found no significant alignments. 

In conclusion, we have found intriguing evidence of non-random alignments in the multipoles £ = 2 — 5, which have 
not only survived, but have indeed been strengthened by the recently released 3-year WMAP data. 
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statistic 


Co-added 1-yr 


Co-added 3-yr 


ILC 1-yr 


ILC 3-yr 


TOHMaske 


TOHkp2 


^22 


0.294 


0.242 


0.242 


0.437 


0.437 


0.294 


S23 


0.011 


0.011 


0.007 


0.021 


0.006 


0.017 


^24 


0.519 


0.637 


0.519 


0.803 


0.774 


0.519 


S25 


0.922 


0.903 


0.937 


0.903 


0.864 


0.913 


S33 


0.457 


0.530 


0.457 


0.549 


0.270 


0.586 


S:u 


0.965 


0.973 


0.978 


0.978 


0.990 


0.965 


S35 


0.957 


0.943 


0.957 


0.921 


0.921 


0.921 


S44 


0.960 


0.960 


0.975 


0.985 


0.905 


0.960 


545 


0.750 


0.750 


0.750 


0.750 


0.750 


0.750 


555 


0.714 


0.714 


0.714 


0.714 


0.664 


0.689 


D23 


0.034 


0.046 


0.031 


0.056 


0.019 


0.051 


D24 


0.771 


0.884 


0.771 


0.922 


0.937 


0.085 


D25 


0.854 


0.786 


0.907 


0.744 


0.550 


0.786 


Dii 


0.268 


0.333 


0.276 


0.351 


0.187 


0.361 


D34 


0.802 


0.888 


0.916 


0.916 


0.965 


0.850 


O35 


0.711 


0.711 


0.797 


0.607 


0.375 


0.607 


D44 


0.913 


0.913 


0.913 


0.913 


0.623 


0.913 


D45 


0.342 


0.342 


0.342 


0.342 


0.786 


0.499 


O55 


0.970 


0.970 


0.989 


0.970 


0.970 


0.970 


R22 


0.828 


0.931 


0.880 


0.647 


0.647 


0.811 


-R23 


0.071 


0.030 


0.091 


0.091 


0.176 


0.115 


-R24 


0.576 


0.767 


0.710 


0.896 


0.817 


0.576 


i?25 


0.823 


0.759 


0.914 


0.759 


0.598 


0.759 


-R33 


0.424 


0.506 


0.424 


0.533 


0.294 


0.533 


-R34 


0.706 


0.902 


0.937 


0.937 


0.986 


0.902 


i?35 


0.834 


0.834 


0.899 


0.834 


0.834 


0.834 


R44 


0.892 


0.939 


0.968 


0.968 


0.485 


0.939 


R45 


0.338 


0.198 


0.198 


0.198 


0.670 


0.338 


R55 


0.977 


0.977 


0.977 


0.977 


0.945 


0.945 


B23 


0.940 


0.920 


0.920 


0.920 


0.968 


0.883 


B24 


0.139 


0.192 


0.065 


0.192 


0.033 


0.117 


B25 


0.344 


0.256 


0.392 


0.256 


0.806 


0.298 


B33 


0.696 


0.582 


0.644 


0.540 


0.749 


0.554 


B34 


0.651 


0.426 


0.344 


0.262 


0.189 


0.506 


B35 


0.591 


0.414 


0.506 


0.414 


0.784 


0.414 


B44 


0.079 


0.116 


0.201 


0.249 


0.300 


0.116 


B45 


0.336 


0.336 


0.497 


0.497 


0.497 


0.336 


B55 


0.444 


0.444 


0.338 


0.338 


0.136 


0.233 



TABLE I: Probabilities P+ for tiie tests S, D, R and B. These are tiie probabilites that a random map would have values for 
the tests S, D, R and B which are higher than the map's values. The probability that a random map has a value lower than 
the map's is simply P_ = 1 — P+. 



Likelihood 


Co-added 1-yr 


Co-added 3-yr 


ILC 1-yr 


ILC 3-yr 


TOHMaskS 


TOHkp2 


L3a{S,D,R,B) 


3.2 X lO"^-'' 


2.0 X 10"" 


7.7 X 10"^** 


3.6 X 10"" 


3.7 X 10"^** 


4.8 X 10"^^ 


L29iS,D,R) 


1.7 X 10"" 


4.3 X 10"^=^ 


5.2 X 10"^^ 


5.1 X 10""' 


3.0 X 10"'^ 


1.3 X 10-"' 


L29iS,R,B) 


5.1 X 10-"' 


4.3 X 10"" 


3.0 X 10"^^ 


8.8 X 10"" 


3.9 X 10"^^ 


3.6 X 10--^ 



TABLE IL Total likelihoods LssiS, D, R, B), L;is{S, D, R) and 1/29(5, R,B). The normalized histograms for the total likelihoods 
found by simulating 25,000 maps are shown in Fig. 6. 
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